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Aims. We report on a multi-band high-energy observing campaign aimed at studying the long term spectral variability of the 
Anomalous X-ray Pulsar (AXP) 1RXS J170849.0-400910, one of the magnetar candidates. 

Methods. We observed 1RXS J170849.0-400910 in Fall 2006 and Spring 2007 simultaneously with SwiftPCRT, in the 0. 1-10 keV en- 
ergy range, and with INTEGRAL/IBIS, in the 20-200 keV energy range. Furthermore, we also reanalyzed, using the latest calibration 
and software, all the publicly available INTEGRAL data since 2002, and the soft X-ray data starting from 1999 taken using BeppoSAX, 
Chandra, XMM-Newton and SwiftPiKT, in order to study the soft and hard X-ray spectral variability of 1RXS J170849.0-400910. 
Results. We find a long-term variability of the hard X-ray flux, extending the hardness-intensity correlation proposed for this source 
over 2 orders of magnitude in energy. 
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CO 1- Introduction 



Anomalous X-ray Pulsars (AXPs, Mereghetti et all |2002|) are 
a peculiar subclass of X-ray pulsars which share the follow- 
ing properties: a narrow range of spin periods (P = 5-12 
s), a typical X-ray luminosity of Lx ~ 10 34 — 10 35 erg s _1 , 
no evidence for Doppler shifts in the light curve, the lack of 
bright optical companions - indicating their isolated nature - 
and a spin down in the range 10~ 13 -10~ 10 s s . Their soft X- 
ray spectra are generally well described by a two component 
model, made of a b lack body with kT ~ 0.4- 0.5 keV, and a 
steep power law (see lWoods & Thomps on 2006, and references 
therein). The detectio n of Soft Gamma-ra y Repeater-like bursts 
from five AXPs (e.g. iGavriil et alj 120021) has pushed forward 
the interpretation of the AXPs as magnetars, i.e. neutron stars 
with an ult ra-high magnetic field (B ~ 10 G). The magne- 
tar model dDuncan & Thompson! 1 19921) . developed initially to 
describe the SGRs phenomenology (bursts, timing, persistent 
emission), seems quite capable to describe the AXPs character- 
istics as well. 

Since the AXPs spectra below ~10 keV are rather soft, the 
first INTEGRAL detections above 20 keV of very hard high- 
energy tails associa t ed with these objects c a me as a surprise 
( Kuiper et al. 2006 1 Iden Hartog et al.l l2006t [Revnivtse v et alJ 
2004; Meregh etti et al.l 120051: iMolkov et al] l2005t iGotz et all 
2006). AXP spectra flatten (T ~ 1) above 20 keV a nd the pulsed 



fraction of some of them reaches up to 100% ( Kuiper et al. 



2006). The discovery of these hard tails provides new constraints 
on the emission models for these objects since their luminosities 
might well be dominated by hard, rather than soft, X-rays. 

1RXS J170849.0-400910 was discovered dur ing the ROSAT 
all sk y survey. The me asure of its period (fSugizaki et al.l 
119971) . period der ivative dKaspi & Gavriill 120031) and general 
X-ray properties dlsrael et alj 19991) . made it an AXP mem - 
ber. Interesting results have been reported by iRea et al.l ([2005 ), 
who claimed a correlation between the soft X-ray flux and 
spectral hardness, with a marginal evidence of the highest and 
har dest spectra being correlated with the AXP glitching activ - 
ity dDairbsso et al.ll2Q03h iKaspi et al.ll200Ct llsrael et al.ll2007l) 
This c o rrelation has been rece ntly confirmed by Campana et ail 
d2007l) . lRea et alj d2007al) . and llsrael etafl d2007l) using further 
Chandra, SwiftPiRT and RXTE data, but always focussing on a 
limited energy interval (i.e. below ~10 keV). 

In this paper, we present our recent multi-w avelength ob- 
serva tion campaign carried out with INTEGRAL dWinkler et al.l 
2003) and Swift dGehrels et al.ll2004l) . These new data, together 
with a re-analysis of all publicly available INTEGRAL, and soft 
X-ray observations of 1RXS J170849.0-400910, allow us to in- 
vestigate the timing and spectral properties of the source over 
a broader energy range. For our timing analysis we made use 
of the solution recently derived by llsrael et al.l (120071) . see Table 
[3] We report a long-term correlation between the soft and hard 
X-ray emission. Results are discussed in the framework of the 
magnetar model. 
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2. Observations and Data Reduction 

2.1. INTEGRAL//B/S observations 

We selected and analyzed all publicly availa ble IBIS 
(INTEGRAL coded mask imager; lUbertinietalfe oCB) pointings 
within 12° from the direction of the source, for a total of 2550 
pointings of 2-3 ks each. In addition, we analyzed our data 
of the Key Programme observation of the Galactic Centre, 
622 pointings performed in Fall 2006, and Spring 2007, since, 
thanks to its large field of view (29°x29°), the source was 
almost always covered during the se observations. Ou r analysis 
is based on data taken with ISGRI dLebrun et al.ll2003l) . the IBIS 
low energy detector array, which is made of CdTe crystals, and 
is working in the 15 keV-1 MeV energy range. 

We processed the data using the Offline Scientific Analysis 
(OSA) software provided by th e INT EGRAL Science Data 
Centre (ISDC, ICourvoisier et al.1 120031) v6.0. We produced 
the images of each pointing in ten energy bands between 20 
and 300 keV. We then added up the images in order to detect 
the source in the 20-70 keV energy range. To investigate the 
broad-band temporal X-ray variability, the observations were 
ganged-up to overlap available soft X-ray data sets. The six data 
segments used in this study are shown in Table Q] From each of 
these segments we derived spectra from the mosaicked images 
in the 10 energy bands. 

Table 1. INTEGRAL observations and results. Errors are at 1 cr 
c.l. 



a P = Public ; KP = Key Programme. 

b In units of counts s -1 on the ISGRI camera. 



2.2. X-ray observations 

We re-analysed all the public available recent data of 
1RXS J170849.0-400910 taken with X-ray imaging telescopes 
in the soft X-rays (i.e. below ~10 keV). The observation times 
and exposures for each data set are reported in Table [2] 

BeppoSAX, XMM-Newton (PN), and Ch andra data were re - 
duced following the procedures described in | lsrael et alj (|2001|) ; 
iRea et al.l (l2003ll2005h : ICampana et all (120071) . but using the lat- 
est available software and calibration files. 

For the Swift/XRT we analysed only data taken in photon 
counting (PC) mode. These were reduced using the FTOOLS 
v6.2 package and CALDB calibration files v20070531 provided 
by the High Energy Astrophysics Science Archive Research 
Center. Standard cleaning and selections were applied and spec- 
tra were extracted using regions of 40" and 80" radius for source 
and background, respectively. Standard response files have been 
used (we used v.9 response matrices). 

We notice that BeppoSAX observed the source also in the 15 - 
300 keV band, with the PDS instrument dFrontera et alJll997l) . 
This non-imaging spectrometer had a field of view of 1.3° 
(FWHM) and the background subtraction was done with a rock- 
ing system, which switched between the source and two off- 
set background regions. The presence of bright contaminating 



sources in both the background pointings prevented us from an- 
alyzing the background subtracted spectra, so we concentrated 
on the timing analysis of these data and extracted barycentered 
events files for the on-source position (see Sect. |3.2| >. 



Table 2. X-ray data observations and results. Errors are at 1 cr 
c.l. 



Instrument 
Name 


OBS date 
MJD 


Exposure 
ks 


1-10 KeV 
Flux" 


Photon 
Index 


SAX LECS/MECS 
SAX LECS/MECS 
Chandra HETG 
XMM-Newton PN 
Chandra CC 
Swift XRT 
Swift XRT 
Swift XRT 


51268 
52138 
52526 
52879 
53189 
53425 
54017 
54182 


26/52 
77/199 
32 
31 
29 
12 
9 
4 


4 40 +U UJ 

-P.-P.3 
4 4*; +0 01 

-0.02 
t q+0.r 
J - -0.2 

3 00 +0 01 

-QOj, 

4 28 +0 " 

^ ,z -°-0.08 

i ae+o.09 
3 26 +012 


2 63 +u - uv 

— Q.Q7 
O Q 8+0.03 

9 QQ+0.12 
Z " J -0.12 

2 77+0.01 
7 74+8:8i 

• ' -Q.Q8 

2 47+0.06 
2 99+8:8§ 
2 76^ : ^ 


a Absorbed flux in units of 10 1 


'erg cirT 2 s 







Obs. 


T Start 


Tj> top 


Exposure 


Type 11 


20-70 keV 


Number 


MID 


MJD 


Ms 




Flux b 


1 


52698.32 


52751.59 


0.83 


P 


0.41+0.06 


2 


52859.63 


52919.88 


1.16 


P 


0.27±0.06 


3 


53226.18 


53298.48 


1.0 


P 


0.42±0.06 


4 


53408.32 


53481.27 


1.5 


P 


0.44±0.05 


5 


53990.59 


54013.77 


0.46 


KP 


0.16±0.09 


6 


54159.36 


54184.07 


0.52 


KP 


0.23±0.09 



3. Analysis and Results 

3. 1 . Spectral analysis 

The spectral parameters reported in Table|2]have been derived by 
fitting all the datasets simultaneously in 1-10 keV energy range 
(except for Chandra data which were limited to 8 keV) by using 
an absorbed black body plus power law model, with a multi- 
plicative constant factor taking into account inter-calibration is- 
sues, not allowed to vary for the same instrument and observ- 
ing mode. While the parameters of the power law have been 
left free to vary, the absorption column density, and the black 
body temperature have been forced to be the same for all in- 
struments. We found a good fit Cf 2 /d.o.f. = 1188/1264 = 0.94), 
and the derived values we re Nh = 1.36(1) x 10 22 cm " 2 (solar 
abundances assumed from lAnders & Grevessd 19891 and pho- 
toelectric cross-section from Balucinska-Church & McCammon 
11992b , and kT = 0.44(1) keV; we verified that, if leaving these 
parameters free to vary, they do not change significantly among 
all the observations. The fluxes and photon indices are reported 
in Tableland Fig.Qffl Our new XRT observation campaign (last 
two points in Fig.[TJ shows that the source entered a new low/soft 
state (similar to the one measured with XMM-Newton in 2003), 
with a flux a factor -1.5 lower than the one measured with XRT 
in 2005. By adding these new points to the long term variabil- 
ity study of t he source, we co nfirm the flux-hardness correlation 
proposed bv lRea et"ai] (|2005j). We also notice that the last three 
points reported in Fig.Q] are particularly compelling since, be- 
ing taken with the same instrument (Swift/XRT), they are not 
affected by cross-calibration uncertainties. In addition, our new 
hard-X data show that the long term variation in flux is correlated 
over more than two orders of magnitude in energy. In fact, IBIS 
observations taken quasi-simultaneously with the last two XRT 
ones (i.e. in late 2006 and early 2007) indicate that the source is 
hardly detected above 20 keV, while the hard X-ray count rates 
measured before followed well the variations measured in the 
soft X-ray range (see Table Q] and Fig.|TJ. Unfortunately, due to 
the faintness of the source we could not statistically prove spec- 
tral changes at high energies, by comparing different INTEGRAL 



1 Note that all the spectral values we report here using the latest cal- 
ibrations and software are consistent with previous findings, except for 
the XMM flux which, using SAS v. 7.1.0, appears slightly higher than 
using SAS 5.4.1 as in Rea et al. (2005); the spectral parameters are in- 
stead consistent. 
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Fig.l. Left plot. Upper Panel: hard X-ray fluxes derived with INTEGRAL/IBIS (20-70 keV). Middle Panel: absorbed 1-10 keV 
fluxes in units of 10~" erg cirT 2 s _1 derived from recent observations of X-ray imaging telescopes as a function of time. Lower 
Panel: photon i ndices measured in the 1-10 keV energy band. Vertical dashed lines mark the times of four observed glitches, (see 
llsrael et alJkOOTl and references therein). Right plot. Folded light curves of 1RXS J170849.0-400910, obtained with RXTEfPCA 
(panels A-D; using the whole 2004 available data) and IBIS data (obs. 3) folded on the RXTE timing solution (panel E, data from 
obs. 3). The energy bands are, from top to bottom: 2.5^1, 4-8, 8-16, 16-32, 20-200 keV. Panels F and G represent the MECS (1-10 
keV) and PDS (20-200 keV) light curves registered during the second (2001) BeppoSAX observation. Note that panels F and G are 
phase-aligned between themselves, but not with the others. Errors are for both plots reported at 1 <x c.l. 



observations. In order to obtain a statistically significant high en- 
ergy spectrum, we co-added the IBIS data from observations 3 
and 4. The resulting 20-200 keV spectrum is well fitted by a sin- 
gle power law, without the need for a cutoff, with photon index 
F = 1.46 + 0.21. The 20-100 keV flux is (3.6 + 0.5) x 10 11 ergs 

-2-1 

cm s . 

In the attempt to characterize the source spectrum in the low 
and high state of the source, we fitted simultaneously the IBIS 
observation 2 with the XMM-NewtonfPN data and the IBIS ob- 
servation 4 with the XRT data. We found that the high energy 
component is well above the extrapolation of the power law 
derived from soft X-ray data (below ~10 keV), while a three 
model component, with an additional power-law describing the 
high energy data points, fits well the entire broadband spectrum 
(1-200 keV). Motivated by this, we tested if the spectral varia- 
tions at low energies could be induced just by the variation of 
the high energy power law. Indeed, by fixing the spectral param- 
eters derived from the XMM-Newton observations, one can fit 
the broad band XMM/INTEGRAL (2003) and Swift/INTEGRAL 
observations (2005) by simply changing the slope and normal- 
isation of the high energy power law. Unfortunately, due to the 
uncertainties in the inter-calibration, and to the low statistic of 
the Swift data, we cannot draw a firm conclusion on this point. 

We also modeled the multi-band spectrum in the "low" 



(2003) and "high" (2005) state with a resonant cyclotron scatter- 
ing model plus a power-law, and we find a good fit in the "low" 
state with B = 0.49(1), kT = 0.36(3), t = 1.2(1) and Y = 1.0(1) 
(see Rea et al. 2007b for details on the model). In the "high" 
state, B and kT do not vary, while we found a r = 1.8(1) and 
F = 0.8(1), which reflects the hardening of the spectrum. The 
"high" state best fit spectrum is shown in Fig. [2] 

3.2. Timing analysis 

IBIS data alone do not allow a blind search for timing signa- 
tures. In order to perform a timing analysis in the hard X-ray 
range, we first derived a phase coherent timing solution using 
publicly available RXTEfPCA data covering the same period 
of the IBIS data. We discovered two new glitches, close to the 
epoch of o ur IBIS observatio ns 3 and 4, which are described 
in detail in llsrael et alj d2007l) . We used the temporal solutions 
derived there (see Table [3]), for the time periods covering the 
INTEGRAL data before and after the first new glitch (and 
prior the second one). We then extracted the IBIS events in the 
20-200 keV band from the pixels that were illuminated by the 
source for at least the 60% of their surface. To optimize the 
signal to noise ratio, we restricted our search to the high state 
periods, namely fall 2004 and spring 2005 (obs. 3 and 4). For 
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Fig. 2. 2005 "high" state phase averaged spectrum of 
1RXSJ170849.0-400910. Upper panel: the points below 10 
keV are taken with Swift/XRT, while above 20 keV IBIS/ISGRI 
data points are shown. The continuous line represents the best fit 
model (RCS + power law) described in the text. Lower panel: 
residuals with respect to the best fit model, plotted as a dashed 
line. 



obs. 3 (obs. 4 results are not reported, but are consistent) we 
folded the data and, by using the Z 2 test we find that the signal 
is pulsed at a ~5.7 cr level (chance probability of 1.6 x 10~ 8 ). If 
we divide the data in two energy bands, 20-60 and 60-200 keV, 
the significance lowers to 3.2 and 4.2 cr, respectively. The pulse 
profile of the source, as a function of energy, is shown in Fig.Q] 
As it can be seen, there is a clear energy dependence of the 
pulse shape morphology: the peak which is predominant in the 
soft band disappears and a secondary peak grows with energy 
becoming the most prominent one above ~8 keV. The pulsed 
fraction cannot be easily derived from the PCA (a non-imaging 
instrument) data due to the uncertainties in the background es- 
timation. Besides, the IBIS data are not easy to handle, because 
of large background variations due the to fact that the source 
is at different off-axis angles during the different pointings. In 
addition, due to the large IBIS field of view, the events can be 
polluted by nearby variable sources. Nevertheless, by averaging 
the background over the entire IBIS observations, we could 
measure a pulsed fraction by dividing the difference between 
the maximum and the minimum of the folded light curves by 
the count rates derived from imaging. The resulting values 
are: -25% (in the 20-60 keV band), -60% (60-200 keV), and 
-40% (20-200 keV). 

We did not find significant differences in the IBIS data 
as far as it concerns the pulsed profiles or the pulsed fraction 
between observations 3 and 4, i.e. before and after the third 
glitch reported in Fig.Q] 



Tabl e 3. T i ming solutions used for obs. 3 and 4, derived in 
llsrael et alj d2007l) 



Obs. 


Epoch 


Period 


Period derivative 


Number 


MJD 


s 


10-" s s' 1 


3 


53189.0 


11.0023066(2) 


1.915(5) 


4 


53189.0 


11.002279(2) 


2.01(1) 



We also reanalyzed the archival PDS data taken with 
BeppoSAX during the pointings listed in Table|2] Using the 
timing solutions derived from th e simultaneous MECS data 
dRea et alj|200l llsrael et all 1200 lb we find, and report for the 
first time, pulsations in the PDS data, in the 20-200 keV energy 
band. Using Z 2 statistics we did not detect the pulsations in the 
first dataset, the shorter one, with high confidence (2.7 cr), while 
in the second dataset they were detected at 5.9 cr. By computing 
the counts in excess of the phase minimum, we found a 20-200 
keV pulsed flux of 0.025 + 0.010 (for the second observation; 
the results for the first one are compatible but less significant). 
By assuming a powe r-law model with photon index T — 1 for 
the pulsed spectrum dKuiper e t al. 2006), this gives a flux of ~ 
0.7 X 10~ n erg cm' 2 s _1 , slightly lower but still compatible with 
that derived from IBIS data by using the same approach, namely 
~ 10" 11 erg cm -2 s _1 . 

4. Discussion and Conclusions 

We have analyzed all the currently available high energy data 
(1-200 keV) of 1RXSJ170849.0-400910 and continued the 
long term variability study of the source in the soft X-ray range 
(^10 keV ) , confi rmin g the flux-hardness corr elation proposed by 
iRea et al.l (120051) and ICampana et all (120071) . Moreover, we re- 
port for the first time the discovery of hard X-ray long term flux 
changes, and show that there is a possible correlation with the 
flux variations detected at lower energies. Thanks to the dense 
monitoring of the last years, we are able to correlate these vari- 
ations with the p resence of the two new glitches (discussed in 
llsrael et aTll2007l) . The only other magnetar for which long term 
variability in this band has been reported is SGR 1806-20, for 
which the spectral hardening app ears to be correlated with a 
burst rate increase dGotz et al.ll2007l). 

As proposed in iThompson et al.l d2002l) . IRea et all d2005l) . 
and ICampana et al.l d2007l) . an increase in the source activity 
(bursts, glitches) and a simultaneous spectral hardening may be 
caused by a growing twist in the magnetosphere; this may be 
also responsible for the transient appearance of a cyclotron line 
during the "high" emission state dRea et al.ll2.Q03l 120051) . If this 
is the case, data presented here support models in which also the 
hard X-ray tails are produced by mechanisms whose strength 
increases with the twist. Indeed, it may also be that the flux vari- 
ations in the hard X-rays dominate and drive those detected at 
lower energies, if the hard X-ray component does not sharply 
cut off below ~10 keV (as it may be the case for a synchrotron 
component). Quite recently, Thompson & Beloborodov (2005) 
discussed how soft gamma-rays may be produced in a twisted 
magnetosphere, proposing two different scenarios: either ther- 
mal bremsstrahlung emission from the surface region heated 
by returning currents, or synchrotron emission from pairs cre- 
ated higher up (~100 km) in the magnetosphere. While both 
scenarios predict a power-law-like spectral distribution for the 
20-100 keV photons, the cut-offs of the high energy emis- 
sion are markedly different in the two cases, 100 keV vs. 1 
MeV. A third scenario involving resonant magnetic Compton 
up-scattering of soft X-ray photons by a non-thermal popu- 
lation of highly relativi stic electrons has been proposed by 
iBaring &~ Harding (2007) . In thi s contest we notice that, as for 
4U 0142+614 dRea et al.ll2007bl) . another member of the AXP 
class, multi-band spectral observations seems to disfavor the 
thermal bremsstrahlung model. In fact, for any choice of tem- 
perature, a model which fits well the data below ~150keV, will 
over-predict by more than a fact or of 5 the C OMPTEL upper 
limits derived in the MeV band (Kuiperet al. 2006). However, 
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these observations are not simultaneous and due to the variabil- 
ity of the source, no firm conclusion can be drawn. 

IBIS data confirm the large pulsed fraction 25%-60%, grow- 
ing with energy, which has been reported earlier at high energies 
for other AXPs (and for 1RXS J1708 49.0-400910using a differ- 
ent data set) bv lKuiper et alj d2006l) . There is a clear indication 
for a phase difference of A<p ~ 0.3 - 0.4 between the maxima of 
the soft (<8 keV) and hard (>8 keV) X-rays light curves mea- 
sured in the 2004 with ISGRI/PCA data (see Fig.Q}. The same 
phase shift was already present in the 2001 MECS/PDS data: 
within the errors, we find no evidence for any significant evolu- 
tion in A<p or pulsed flux, despite the three glitches that occurred 
between the observations . The phase shift , similar to the one de- 
tected at lower energies dRea et al.ll2.Q03h , still lacks a clear in- 
terpretation in the magnetar framework, and it may indicate that 
the location of production of the two (soft/hard) component is 
different, but stable with time. 
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